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Abstract 

In this paper, we develop several related finite dimensional variational principles for discrete optimal 
transport (DOT), Minkowski type problems for convex polytopes and discrete Monge-Ampere equation 
(DMAE). A link between the discrete optimal transport, discrete Monge-Ampere equation and the power 
diagram in computational geometry is established. 

1 Introduction 

1.1 Statement of results 

The classical Minkowski problem for convex body has influenced the development of convex geometry and 
differential geometry through out the twentieth century. In its simplest form, it states, 




Figure 1 : Minkowski problem 

Problem 1 (Minkowski problem for compact polytopes in R n ) Suppose ni, n& are unit vectors which 
span R n and A\, > so that Y2i=i ^i n i = 0- Find a compact convex polytope P C R n with exactly 
k codimension-1 faces F\, Fk so that rij is normal to F$ and the area of Fi is A{. 

Minkowski's famous solution to the problem says that the polytope P exists and is unique up to parallel 
translation. Furthermore, Minkowski's proof is variational and suggests an algorithm to find the polytope. 

Minkowski problem for unbounded convex polytopes was solved by Alexandrov in his influential book 
on convex polyhedra [T). In particular, he proved the following fundamental theorem (Theorem 7.3.2) which 
is one of the main focus of our investigation. 

Theorem 1.1 (Alexandrov ) Suppose f£ is a compact convex polytope with non-empty interior in R n , pi , . . . , p\. C 
R n are distinct k points and A\,...,A^ > so that Y2i=i Ai = woZ(O). Then there exists a vector 
h = (hi, hk) € R k , unique up to adding the constant (c, c, ...,c), so that the piecewise linear con- 
vex function 



u(x) = max{x • pi + hi} 



satisfies vol({x € £l\Vu(x) = pi}) = Ai 
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Figure 2: Discrete Optimal Transport Mapping (left to right): map W% to pj. Discrete Monge-Ampere 
equation (right to left): vol(Wi) is the discrete Hessian determinant of Pi. 

The functions u and \/u(x) in the theorem will be called the Alexandrov potential and Alexandrov map. 
Alexandrov's proof is non-variational and non-constructive. Producing a variational proof of it was clearly 
in his mind. Indeed, on page 321 of 1> he asked if one can find a variational proof and considered such 
proof "is of prime importance by itself. One of the main results of the paper gives a (finite dimensional) 
variational proof Alexandrov's Theorem 11.11 Indeed, we give a variational proof of a general version of 
Theorem 11.11 (Theorem 11.21 below) and produce an algorithm for finding the function u. In recent surge 
of study on optimal transport, Theorem 11.11 is reproved and is a very special case of the seminal work 
of Brenier (see for instance [9|, Theorem 2.12(h), and Theorem 2.32). Brenier proved that the function 
yu minimizes the quadratic cost J n \x — T(x)\ 2 dx among all measure preserving maps (transport maps) 
T : (Q, dx) — > (R n , Yli=i Ai5 Pi ). Here 5 P is the Dirac measure supported at the point p. Thus our work 
produces a variational principle and an algorithm for finding Alexandrov maps with finite images. 

1.2 Variational principles 

Here is a simple framework which we will use to establish variational principles for solving equations in this 
paper. Suppose X C R k is a simply connected open set and A(x) = (Ai(x), ...,Ak(x)) : X — > R k is a 
smooth function so that ^gff = dA Q^ for all Then for any given B = (B\, B^) £ R k , solutions x 
of the equation A(x) = B are exactly the critical points of the function E(x) = J x Yli=i(Ai(x) ~ Bijdxi. 
Indeed, the assumption d ^ x ^ = ^gj^ says the differential 1-form uj = J2i=i(Ai(x) — Bi)dxi is closed 
in the simply connected domain X. Therefore the integral E(x) = ff u: is well defined independently of 

the choice of the path from a to x. By definition, = Ai{x) — B{, i.e., xjE(x) = A(x) — B. Thus 

A{x) = B is the same as \/E{x) = 0. 

All variational principles established in this paper use the above framework. We will use the above 
framework to give a variational proof of Alexandrov's theorem. 

The paper will mainly deal with piecewise linear (PL) convex functions. Here are the notations. Given 
pi, ...,Pk € R n and h = (hi, h^) € R k , we use u(x) = Uh{x) to denote the PL convex function 

u h [x) = max{x -pi + hi}, 

i 

where u ■ v is the dot product. Let Wi(h) = {x S R n | V u i x ) = Pi} be the closed convex polytope. It is 
well known that Wi(h) may be empty or unbounded. One of the main result we will prove is, 
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Theorem 1.2 Let Cl be a compact convex domain in R n and {pi, ...,pk} a set of distinct points in R n and 
a : Cl — > R be a positive continuous function. Then for any A\, > with Yli=i = In o~(x)dx, 
there exists b = (b\, bjS) G R k , unique up to adding a constant (c, c), so that J w .^ n Q o~(x)dx = 
for all i. The vectors b are exactly maximum points of the convex function 

ph k» k 

E(h)= / V / a(x)dxdhi -S^hiAi (1) 

J a i=1 Jw l (h)nn i=l 

on the open convex set H = {h € R k \vol(Wi(h) DCl) > Ofor all i}. In fact, E(h) restricted to Hq = 
HC\{h\ X^i=i hi = 0} is strictly convex. Furthermore, V 'u\, minimizes the quadratic cost f n \x — T(x)\ 2 adx 
among all transport maps T : (CI, adx) — > (R n , Yli=i Ai§ Pi ). 

We remark that Alexandrov's theorem corresponds to a = 1. The existence and the uniqueness of 
Theorem 11.21 are special case of important work of Brenier on optimal transport. Our main contribution is 
the variational formulation. The Hessian of the function E(h) has a clear geometric meaning and is easy 
to compute (see equation [6]), which enables one to efficiently compute Alexandrov map using Newton's 
method. Furthermore, as a consequence of our proof, we obtain a new proof of the infinitesimal rigidity 
theorem of Alexandrov that VE : H -> W = {(A ± , A k ) G R k |Aj > 0, ^=1 A i = In o-(x)dx} is a 
local diffeomorphism (see Corollary 13.2) . We remark that Aurenhammer et al. |0 1 also noticed the convexity 
of the function E, and they gave an elegant and simple proof of v u 6 minimizing quadratic cost. 



1.3 Discrete Monge-Ampere equation (DMAE) 

Closely related to the optimal transport problem is the Monge-Ampere equation (MAE). Let Cl be a compact 
domain in R n , g : dCl — » R and A:l!xRx R n — > R be given. Then the Dirichlet problem for MAE is 
to find a function w : Cl — > R so that 

det(Hess(w))(x) = A(x, w(x), yw(x)) 
w\dn = g 

There are vast literature and deep results known on the existence, uniqueness and regularity of the 
solution of MAE. We are interested in solving the discrete version of MAE in the simplest setting where 
A(x, w, Sjw) = A(x) : Cl — > R so that A(Cl) is a finite set. By taking Fenchel-Lengendre dual of the 
Alexandrov potential function u, we produce a finite dimensional variational principle for solving a discrete 
Monge-Ampere equation. 

In the discrete setting, one of the main tasks is to define the discrete Hessian determinant for piecewise 
linear function. We define, 

Definition 1.3 Suppose (X, T) is a domain in R n with a convex cell decomposition T and w : X —■ ~Ris a 
convex function which is linear on each cell (a PL convex function). Then the discrete Hessian determinant 
of w assigns each vertex vofT the volume of the convex hull of the gradients of the w at top-dimensional 
cells adjacent to v. 

One can define the discrete Hessian determinant of any piecewise linear function by using the signed vol- 
umes. This will not be discussed here. With the above definition of discrete Hessian determinant, following 
Pogorelov [7], one formulates the Dirichlet problem for discrete MAE (DMAE) as follows. 

Problem 2 (Dirichlet problem for discrete MAE (DMAE)) Suppose Cl = conv(v\, v m ) is an n- 
dimensional compact convex polytope in R n so that Vi ^ conv(v\, Vi-\, fj+i, Vk) for all i. Let 
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pi,...,Pk beinint(Q). Given any g\, g m € R and A±, Ak > 0, find a convex subdivision T of Q with 
vertices exactly {vi, v m ,p\, ...,Pk} and a PL convex function w : f2 — > R which is linear on each cell of 
T so that 

(a) (Discrete Monge-Ampere Equation) the discrete Hessian determinant of w at pi is A{, 

(b) (Dirichlet condition) w(vi) = gi. 

In Q, Pogorelov solved the above problem affirmatively. He showed that the PL function w exists and 
is unique. However, his proof is non- variational. We improve Pogorelov 's theorem to the following. 

Theorem 1.4 Suppose $7 = conv(v\, ...,v m ) is an n-dimensional compact convex polytope in R n so that 
Vi $l conv(v±, v Vi+i, Vk) for all i and p\, ...,pk are in the interior of Cl. For any g\, g% € R 
and A±, A}. > 0, there exists convex cell decomposition T having and pj as vertices and a piecewise 
linear convex function w : (£1, T) — > R so that w(vi) = gi,i = 1, m and the discrete Hessian determi- 
nant of w at pj is A j, j = 1, Infact, the solution w is the Legendre dual of max{x-pj+hj,x-Vi—gi\j = 
1, k,i = 1, m} and h is the unique critical point of a strictly convex function. 

The paper is organized as follows. In §2, we recall briefly some basic properties of piecewise linear 
convex functions, their dual and power diagrams. Theorems II .2l and ll .4l are proved in §3 and §4. 

2 Preliminary on PL convex functions, their duals and power diagrams 

We collect some well known facts about PL convex functions, their Legendre-Fenchel duals and their rela- 
tions to power diagrams in this sections. Most of the proofs are omitted. See Aurenhammer [2 1, Passera and 
Rullgard |6], Siersmas and van Manen [8 ] and others for details. 

The following notations will be used. For u, v,p%, ...,Pk £ R n , we use u ■ v to denote the dot product 
of u,v and conv(p\, ...,Pk) to denote the convex hull of {p±, ...,Pk} C R n . A convex polyhedron is the 
intersection of finitely many closed half spaces. A convex polytope is the convex hull of a finite set. The 
relative interior of a compact convex set X will be denoted by int(X). 

2.1 Legendre-Fenchel dual and PL convex functions 

The domain of a function / : R n — >(— oo,oo], denoted by D(f), is the set {x S R n |/(x) < oo}. A 
function / is called proper if D(f ) ^ 0. For a proper function / : R n — > (— oo, oo], the Legendre-Fenchel 
duality (or simply the dual) of / is the proper function /* : R n — >(— oo, oo] defined by 

r(y) = sup{x-y-/(x)|xeR n }. 

It is well known that /* is a proper, lower semi continuous convex function. For instance, for the linear 
function f(x) = a ■ x + b, its dual /* has domain D(f*) = {a} so that /*(a) = —b. The Legendre-Fenchel 
duality theorem says that for a proper lower semi continuous convex function /, (/*)* = /. 

For P = {pi, ...,pk} C R n and h = (hi,...,hk) £ R k , we define the piecewise linear (PL) convex 
function Uh (x) to be 

u(x) = (u h (x) = Uh,p{x)) = maxjpj • x + hi\i = 1, k} (3) 
The domain D(u*) of the dual u* is the convex hull conv(pi, ...,Pk) so that 

k k k 

u*(y) = min{~y^ tjh^tj > 0, ^ U = 1,^^ = y} (4) 

i=l i=l t=i 

(See theorem 2.2.7 of Hormander's book on Notions of Convexity [5 ]). In particular, u* is PL convex in the 
domain D(u*). For instance if h = 0, then u*(y) = on D(u*). Another useful consequence is, 
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Figure 3: PL-convex function and its induced convex subdivision. 



Corollary 2.1 If pi £ conv(pi, ...,Pi-x,p i+1 , ...,Pk), then 



u*{pi) = -hi. 



(5) 



Indeed, the only way to express pi as a convex combination of p\, ...,Pk is p% = 1 • p%- Thus © holds. 
2.2 PL convex functions, convex subdivisions and power diagrams 

A PL convex function / defined on a closed convex polyhedron K produces a convex subdivision (called 
natural subdivision) T of K. It is the same as the power diagram used in computational geometry. Let us 
recall briefly the definition (see for instance [6]). A convex subdivision of K is a collection T of convex 
polyhedra (called cells) so that (a) K = U^g-pr, (b) if a, r G T, then a n r G T, and (c) if a G T and 
r C a, then r G T if and only if r is a face of a. The collection T is determined by its top-dimensional 
cells. The set of all zero-dimensional cells in T, denoted by T°, is called the vertices of T. 

If / is a PL convex function defined on a convex polyhedron K, the natural convex subdivision T of K 
associated to / is the subdivision whose top-dimensional cells in T are the largest convex subsets on which 
/ are linear. The vertices of / are defined to be the vertices of T. Suppose {vi, ...,v m } is the set of all 
vertices of /. Then / is determined by its vertices {v{\ and the values at the vertices {f(v{)}. Indeed the 
graph of / over K is the lower boundary of the convex hull conv({v\, f(vi)),..., («/(%))). Recall that if 
P is a convex polyhedron in R n x R, then the lower faces of P are those faces F of P so that if x G F, 
then x — (0, 0, A) is not in P for all A > 0. The lower boundary of P is the union of all lower faces of P. 

One can also describe T by using the epigraph. The epigraph {(x,t) G K x R|i > f(x)} of / is 
naturally a convex polyhedron. Each cell in T is the vertical projection of a lower face of the epigraph. 

Since the dual function /* is also PL convex on its domain D(f*), there is the associated convex sub- 
division T* of D(f*). These two subdivisions (D(/),T) and (D(f*),T*) are dual to each other in the 
sense that there exists a bijective map T — > T* denoted by a — >• a* so that (a) a, r G T with r C a if 
and only if a* C r* and (b) if r C a in T, then the cone(r,a) is dual to cone(a* , r*). Here the cone 
cone(r, a) = {t(x — y)\x G <r, y G r, i > 0} and dual of a cone C is {x G R n | y • x < for all y G C}. 
See proposition 1 in section 2 of Passera and Rullgard [6 ]. 

For the PL convex function Uh{x) given by ©, define the convex polyhedron W{ = Wi(h) = {x G R n [ 
x ■ pi + hi > x ■ pj + hj for all _?'}. (Note that Wi may be the empty set.) By definition, the con- 
vex subdivision T of R n associated to Uh is the union of all Wi's and their faces. Identity (0]) for u* 
says that the graph {(y,u*(y))\y G conu(pi, ...,£>&)} of u* is the lower boundary of the convex hull 
conv((pi, -hi), (p fc , -/j-j))- 
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We summarize the convex subdivisions associated to PL convex functions as follows, 

Proposition 2.2 (a)Ifint(Wi(h)) ^ and p\, ...,pk are distinct, then int(Wi(h)) = {x G H n \x-pi + hi > 

maxj^i{x ■ pj + hj}}. 

(b) If pi ^ conv(pi, ...,pi-\,pi + i, —,Pk), then int(Wi(h)) ^ and Wi(h) is unbounded. 

(c) If conv(pi, ...,pi_i,pi + i, ...,Pk) is n-dimensional and pi G int(conv(pi, ...,pk)), 
then Wi(h) is either bounded or empty. 

(d) Ifpi, ■■■,Pk are distinct so that int(Wi(h)) ^ for all i, then the top-dimensional cells of T associ- 
ated to Uh are exactly {Wi(h)\i = 1, k}, vertices of u* and the dual subdivision T* of conv(p\, —,Pk) 
are exactly {p±, ...,pk}. 

(e) For any distinct pi, ...,pk G R n , there is h G R k so that int(Wi(h)) ^ %for all i. 

Proof To see (a), by definition {x G R n |x • pi + hi > maxj^jjx • pj + hj}} is open and is in Wi(h). 
Hence it is included in the interior int(Wi{h)). Let Lj(x) = x ■ pj + hj. By definition, Li(x) > Lj(x) 
for all x G Wi(h). It remains to show for each p G int(Wi), Li(p) > Lj(p) for j ^ i. Take a point 
q G int(Wi(h)) so that Li(q) > Lj(q) (this is possible since Lj ^ Lj for j ^ i). Choose a line segment / 
from q to r in int{Wi(h)) so that p G int(I). Then using Li(q) > Lj(q), Li(r) > Lj(r), p = tq + (1 — t)r 
for some t G (0, 1) and linearity, we see that Lj(p) > Lj(p). This establishes (a). 

To see part (b) that int(Wi) / 0, by the identity © for u* hpi with P' = {p\,—,Pk} - {Pi}, we 
have u* hP i{pi) = oo. But u* h P ,(pi) = sup{x • pi — Uh,p'{x)\x G R n }. Hence there exists x so that 
x ■ pi + hi > Uh,P'(x) = maxj^i(x ■ pj + hj), i.e., int(Wi(h)) / 0. Furthermore, u* h p/ (pi) = oo implies 
Wi(h) is non-compact, i.e., unbounded. 

To see (c), suppose otherwise, then the set Wi(h) contains aray {tv + a\t > 0} for some non-zero vector 
v. Therefore, (tv + a) -pi + hi > (tv + a) ■ pj + hj for all j ^ i. Divide the inequality by t and let t — > oo, 
we obtain v ■ pi > v ■ pj for all j ^ i. This shows that the projection of pi to the line {tv\t G R} is not in 
the interior of the convex hull of the projections of {p\, ...,Pk} — {Pi}- This contradicts the assumption that 
Pi is in the interior of the n-dimensional convex hull. 

The first part of (d) follows from the definition. The duality theorem (proposition 1 in section 2 of ||6]) 
shows the second part. 

To see part (e), let us relabel the set pi, ...,Pk so that for all i if j > i then pj is not in the convex hull of 
{pi, ...,pi}. This is always possible due to the assumption that pi,...,Pk are distinct. Indeed, choose a line 
L so that the orthogonal projection of p/s to L are distinct. Now relabel these points according to the linear 
order of the projections to L. 

For this choice of ordering of pi , . . . , p^, we construct h\,... : hj t inductively so that Wi (h) contains a non- 
empty open set. Let h\ = 0, since p2 ^ p\, for any choice of fo, both vol({x\ V u (hi,/t 2 )( x ) = Pi}) > f° r 
i = 1,2. Inductively, suppose hi, hi have been constructed so that v ol({x\ V u (hi,...,hi)( x ) = Pj}) > 
for all j = 1, 2, To construct /ij+i, first note that since pi + \ is not in the convex hull of p\, ...,pi, by 
part (a), for any choice of /ij+i, vol(Wi + i{h\, hi + \)) > and Wi + i(h\, /ij+i) is unbounded. Now 
by choosing hi+i very negative, we can make all vol{Wj{h\, > for all j = 1, 2, + 1. 

It is known that convex subdivisions associated to a PL convex function Uh(x) on R n are exactly the 
same as the power diagrams. See for instance G], 10. We recall briefly the power diagrams. Suppose P = 
{pi , . . . , pk } is a set of k points in R n and w\ , . . . , are k real numbers. The power diagram for the weighted 
points {(pi, W\), (pk,Wk)} is the convex subdivision T defined as follows. The top-dimensional cells 
are U{ = {x G R n ||x — pi\ 2 + Wi < \x — pj\ 2 + Wj for all j}. Here \x\ 2 = x ■ x is the square of the 
Euclidean norm and \x — Pi\ 2 + Wi is the power distance from x to (pi,Wi). If all weights are zero, then 
T is the Delaunay decomposition associated to P. Since \x — pi\ 2 + Wi < \x — pj\ 2 + Wj is the same as 
x ■ x — 2x ■ pi + \pi\ 2 + Wi < x ■ x — 2x ■ pj + \pj\ 2 + Wj which is the same as x ■ Pi — ^(\pi\ 2 + Wi) > 
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Figure 4: Power diagram and its dual weighted Delaunay triangulation. 



x-pj — \{\pj\ 2 +Wj). WeseethatJTj = {x € H n \x-pi+hi > x-pj + hj for all j} where hi = —\{\pi\ 2 +Wi)- 
This shows the well known fact that, 

Proposition 2.3 The power diagram associated to {(pi,Wi)\i = 1, k} is the convex subdivision associ- 
ated to the PL convex function defined by (0) where hi = —\eA-+2£±_ 

2.3 Variation of the volume of top-dimensional cells 

The following is the key technical proposition for us to establish variational principles. 

Proposition 2.4 Suppose a : O — > R is continuous defined on a compact convex domain f] C R n . 
If pi, ...,Pk € R n are distinct and h € R k so that vol(Wi{h) CI O) > for all i, then Wi(h) = 
fw(h)nQ cr ( x )d' x is ci differentiable function in h so that for j ^ i and Wi(h) D and Wj(h) D $7 share a 
codimension-1 face F, 



Proof The proof is based on the following simple lemma. 

Lemma 2.5 Suppose X is a compact domain in R n , / : X — > R is a non-negative continuous function and 
r(x,t) : {(x,t) £ X x R|0 < t < f(x)} — > R is continuous. For each t > 0, let f t (x) = min{t, f(x)}. 
Then W(t) = Jxilo*^ r ( x ' s )ds)dx satisfies 




(6) 



where dA is the area form on F and partial derivative is zero otherwise. In particular, 



dwi(h) dwj(h) 



dhj dhi 




(7) 



and 




(8) 



It is differentiable at to if and only if J, 



{xex\f{x)=t } 



t{x, to)dx = 0. 



The conditions in the lemma can be relaxed. 
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Figure 5 : Variation of the volume of top-dimensional cells 

Proof LetG t (x) = f Q ft{x) t(x, s)ds and M be an upper bound of \t(x, t) j in its domain. Since | min(a, b) - 
min(a, c)| < \b — c|, we have \ft{x) — ft'(x)\ < \t — t'\. Now, for any t / t' , 

.G t (x)-G t ,(x). 1 . [f*& . ^ . M 



t-t' ' \t-t' 



fJt(x, M 

/ a(x,s)ds\ < r—jrlMx) - f t ,(x)\ < M. (9) 



Fix to an d x € X. If f(x) < to, then for t very close to to, Gt(x) = Jq t(x, s)ds. Hence 
lim 4 ^ Gt ^_^o^ x ) = o. If f[x) > to, then for t very close to to, Gt(x) = J* * r(x, s)ds. Hence 
lim t ^t Gt — t — = limf^ J t * T(x,s)(is = r(x,to)- If /0&) = *o> then the above calculations 

shows lim t _ >t + Gt ^_^ to — = and lim t ^ t - = r(x, to)- Therefore, by Lebesgue dominated 

convergence theorem, we have 

am EfiLiZM _ lim / <fMz3M ix _ / T(I , to)& (10) 



and 



lim ffw-ffw _ lim r cw-ci^ _ /■ t(Ij 

i->*0 * — *0 t-+*o 7X * — *0 ■'{as|/(a:)>to} 



This establishes the lemma. 

Fix a < b, we call the domain {(x, t) £ X x R|a < /(x), a < t < mm(/(x), b)} a cap domain with 
base {x|/(x) > a} and top {x\f(x) > b} of height (b — a) associated to the function /. 

To prove the proposition 12.41 let hi = (hi, hi — 5, /ij+i, hj~). For small positive 5 > 0, by 
definition, Wi(h') c W t (h) mdWj(h) c Wj(h'). lfWi(h)nWj(h)nn = 0, then Wj(h)nSl = Wj(h')nn 

for small 5. Hence dV ^jj h ^ = 0. If W{ (h) n SI and VFj (ft)nfi share a codimension- 1 face F, then the closure 
cl(Wj(h') — Wj(h)) is a cap domain with base F associated to a convex function / defined on F. The 
height of the cap domain is | p .^ p .| ^ an d / is PL convex so that the (n — 1) -dimensional Lebesgue measure 
of set of the form {x € F\f(x) = t} is zero. Furthermore for 5 > 0, by definition 



Wj(h') — Wj(h) 1 
5 = 5 



a(x)dx 
Wj{h')n(i-Wj(h)nn 



I r rMv) 

$ J Jo T ^' s ) dsd y 
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where y G F is the Euclidean coordinate and r(y, s) is cr expressed in the new coordinate. Thus, by lemma 
12.51 we see 

Hm ^m^m = [ a \ FdA . (12) 

The same calculation shows that for 5 < and close to 0, using the fact that cl(Wj(h) n U — Wj(h') n O) 
is cap with top F, we see that (fT2l holds as well. Finally, if Wi(h) n SI and Wj(/i) n SI share a face of 
dimension at most n — 2, then the same calculation still works where the associate cap domain has either 
zero top area or zero bottom area. Thus the result holds. 



3 A proof of Theorem IL2 



Our proof is divided into several steps. In the first step, we show that the set H = {h G R k |i> ol(Wi(h) n 
SI) > for all i} is a non-empty open convex set. In the second step, we show that E{h) = f n Uh(x)a(x)dx 
is a C 1 -smooth convex function on H so that = Jw (h)nQ a ( x )dx- In the third step, we show that 

E(h) is strictly convex in Hq = H n {h\ Yli=i hi = 0}. In the fourth step, we show that the gradient map 
XjE : H -4 W = {(Ai, ...,A k )\Ai > 0,X)*=iA = Jn a ( x ) dx } is a diffeomorphism. Finally, for the 
completeness, we include a simple proof by Aurenhammer et al. (Lemma 1 in O) to show that Vub is an 
optimal transport map minimizing the quadratic cost. 



3.1 Convexity of the domain H 

We begin with a simple observation that a compact convex set X C R n has positive volume if and only if 
X contains a non-empty open set, i.e., X is n-dimensional. Therefore, vol(Wi(h) D f2) > is the same as 
Wi{h) fl £1 contains a non-empty open set in R n . The last condition, by the above proposition 12.21 a), is the 
same as there exists x 6 f2 so that x ■ Pi > maxj^jjx • pj +hj}. 

Now to see that H is convex, since H = n% =l Hi where Hi = {h G Y&\vol{Wi{h) nH)> 0}, it 
suffices to show that H L is convex for each i. To this end, take a, f3 £ Hi and t G (0, 1). Then there exist 
two vectors v\ , i>2 G O so that t> i -pj + a» > u i -pj + ay and v^-Pi+^i > V2 -Pj + /9j for all j 7^ i Therefore, 
(tvi + (1 - t)v 2 ) -Pi + (tai + (1 - OA) > + (1 - t)v 2 ) -Pj + {t aj + (1 - t)/3j) for all j / i. This shows 
that ta + (1 — is in ifj. Furthermore, each ifj is non-empty. Indeed, given hi, /ij+i, hk, 
by taking /ij very large, we see that h = (hi, %) is in H t . Also, from the definition, Hi is an open set. 
Therefore, to show H is an open convex set, it remains to show that H is non-empty. 

To see H ^ 0, it suffices to show that there exists h so that vol(Wi(h)) > (which could be 00) 
for all i. Indeed, after some translation, we may assume that is in the interior of Q. Due to the fact 
that u Xh (x) = Xu h {x/\) and Wi(Xh) = Wi(h)/X, if vol(Wi(h)) > for all i, then for A > large, 
vol(Wi(\h) n Q) > 0. Now by proposition I2.2f e). we can find h so that vol(Wi(h)) > for all i. Thus 
H + 0. 



3.2 Convexity of the function E{h) and its gradient 

We show that E(h) = f n Uh(x)a(x)dx is convex in h and satisfies dE(h)/dhi = /^ (fe)nO a ( x )dx. Since 
u/j(a;) is the maximum of a collection of functions x ■ pi + hi which are linear in (x,h) G R n x R k , 
Uh(x) is convex in (x, h) G R n x R k . Thus by the assumption that a(x) > 0, Uh{x)a{x)dx is convex 
in h. For a convex function /, the directional derivative of / at the point a in direction d, denoted by 
f'(a,d) = lim t _ >0 + ^ a+td ^^ a ^ , always exists in (0, 00]. Furthermore, it is known (proposition 2.3.2 in 
HI ) that if fi are convex and f(x) = maxj =1 {/j(a;)}, then f'(a,d) = m.aXi£K a {fi(<x, d)} where K a = 
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{j\fj(a) = f(a)}. In our case, x ■ Pi + hi is smooth so that d(a>pH-M = Therefore, we obtain 

du h (x) 

——— = max {dij} 
otij ieK xJl 

where K x ^ = {m\x • p m + h m = Uh(x)}. By proposition \2.2\ a), if x G int(Wj(h)), then 9u q^ = 1 and 

if x € mt(Wi(/i)) for % ^ j, then dU g^ = 0. In particular, this shows \uh(x) — Uh>(x)\ < \h — h'\. 
By the dominant convergence theorem in real analysis and continuity of a(x), we obtain 

dE(h) f du h (x) 



dhi 



f — jf- — a(x)dx = f a(x)dx 



This shows that E(h) = J h Yli Iw (h)nn (T (x)dxdhi + c for some constant c. 
3.3 The strict convexity of E(h) 

We will show that the Hessian matrix of E(h) has a 1 -dimensional null space spanned by the vector 

(h)nn a { x )dx. By the calculation above, dE(h)/ dhi — Wi(h) and Yli=i w i(h) 
j Q a(x)da. By proposition 12.41 for i ^ j, dV ^\^ < 0. Furthermore, if Wi(h) and Wj(h) share a 

codimension- 1 face F in S7, then 9V ^^ = — \ p ._ p .\ Jp &\F(x)dA < where dA is the area form on 
the codimension- 1 face F. This implies, 

Corollary 3.1 The Hessian matrix Hess(E) ofE(h) is positive semi-definite with 1 -dimensional null space 
generated by (1, 1, .., 1). In particular, E\h q '■ Hq = {h G H\ Y^i=i hi = 0} — > R is strictly convex. 

Proof By proposition 12.41 and the fact that Yli=i^ w i/^hj = d(vol(Q))/dhj = 0, it follows that the 
Hessian matrix Hess(E) = [ d ^ h ] = [dij] is diagonally dominated (i.e., all diagonal entries are positive 
and all off diagonal entries are non-positive so that the sum of entries of each row is zero). Therefore, it is 
positive semi-definite with kernel containing the vector (1, 1, 1). To see that the kernel is 1 -dimensional, 
supposes = (i>i,..., Vk) is anon-zero vector so that Hess(E)v t = (v* is the transpose of v). Let us assume 
without loss of generality that jttjj = maxj{|t;j|} and Vi x > 0. In this case, using a^v^ = Ylj^ a hj v j 
and a^jj = — Y^j^tix \ a hj\> we see tnat v j = v h f° r an indices i with a^j / 0. It follows that index set 
I = {i\vi = maxj{\vj\}} has the following property. If i\ G / and ^ 0, then 12 G /. We claim 
that / = {1, 2, k}, i.e., v = ui(l, 1, 1). Indeed, for any two indices i ^ j, since Q is connected 
and W r (h) are convex, there exists a sequence of indices i\ = 1,12, ■■■,i m = j so that Wi s (h) H Q and 
Wi 3+1 (h) H f2 share a codimension-1 face for each s. Therefore ai s i 3+1 / 0. Translating this to the Hessian 
matrix Hess(E) = [ay], it says that for any two diagonal entries an and ajj, there exists a sequence of 
indices i\ = i, 12, i m = j so that ai s : ls+1 < 0. This last condition together with the property of I imply 
I = {1, 2, k}, i.e., dim(Ker(Hess(E))) = 0. 



3.4 VE is a local diffeomorphism 

With above preparations, we can show the gradient map VE is a diffeomorphism. Let $ = VE : Hq — > 
TV = {(At, G R k |A > 0,X]f = iA = Jn^i^dx} be the map sending h to {wi(h), w k {h)). 
By the calculation above, $(/t) is the gradient of E\h at h. Since ^|h has positive definite Hessian matrix 
in Hq, its gradient $ is an injective local diffeomorphism from Hq to W. But dim{Ho) = dim{W), thus 
&(Hq) is open in Thus to finish the proof that $(i?o) = using the fact that W is connected, we only 
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need to show that <&(Hq) is closed in W. To see this, take a sequence of point M m ) in Hq so that <£(/i( m )) 
converges to a point a G W. We claim that a G $(iTo))- First note that h^'s are bounded in R k . Indeed, 
if not, we can choose a convergent subsequence, still denoted by M" 1 ) , so that there are two indices i ^ j 
with — > oo and hj — oo. Since $7 is compact, we see for m large x ■ pi + /ij"^ > x • + hj 
for all 3; G SI. This shows that Wj{hS m ^) flfl = for m large which contradicts the assumption that 
vol{Wj{h^) nfi) > for m large. As a consequence, we can choose a convergence subsequence, still 
denoted by — > h G R k . For this /i, by the continuity of the map h — > (w\(h), ...,Wk(h) on R k , we 
have = q, i.e., h G H and a G ^(i^o)- 

As a consequence of the proof, we also obtained a new proof of the infinitesimal rigidity theorem of 
Alexandrov. 

Corollary 3.2 (Alexandrov) The map \jE : Hq — > W sending the normalized heights h to the area vector 
(wi(h), Wk(h)) is a local diffeomorphism. 

3.5 V«j is an optimal transport map 

We reproduce the elegant proof by Aurenhammer et al [3 ] here for completeness. Notice the quadratic 
transport cost of Vu^ is J2i=i Iw — x\ 2 a(x)dx. From Proposition 12.31 {Wi,...,Wfc} is the power 
diagram associated to {(pi,Wi = —\pi\ 2 — 26j)} in R d . Suppose {Ui, Uk} is any partition of R d so that 
J v a(x)dx = J w a(x)dx, i = 1, 2, • • • , k. By the definition of the power diagram, we have 

y~] j (\x - pi\ 2 + Wi)a(x)dx <y~] (\x - pi\ 2 + Wi)a(x)dx. 
l=1 JWi i=1 JVi 



Thus 



k n k n 

y2 \x~ Pi\ 2 o-(x)dx) < V / 
i=l J Wi i=1 JUi 



\x — pi\ 2 a(x)dx. 



This shows that Vu;, minimizes the quadratic transport cost. 



4 A proof of Theorem IL4 



We fix #1, —,g m through out the proof. For simplicity, let Pk+j = Vj and h^+j = —gj for j = 1, m and 
let 

Wi (h) = {x G R n [x ■ pi + hi > x ■ pj + hj,j = 1, k + . 

Define 

H = {he Y&\vol{Wi{h)) > 0,* = 1, + m}. 

Lemma 4.1 (a) His a non-empty open convex set in R k . 

fftj For eac/j h £ H and i = l,...,k and j = l,...,m, Wi(h) is a non-empty bounded convex set and 
Wk+j{h) is a non-empty unbounded set. 

Proof The proof of convexity of H is exactly the same as that of §2.1. We omit the details. Also, by 
definition H is open. To show that H is non-empty, using proposition |2.2f e), there exists h G R k so that for 
all i = 1, k, vol(Wi(h)) > 0. We claim for t > large the vector h = h + (t, t) G H. Indeed, let B 
be a large compact ball so that B H Wi(h) 7^ for all i = 1, A;. Now choose t large so that 

min{x • pi + hi\i = 1, 2, k} > max{x • Vj + gj\j = 1, m}. 
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For this choice of h, by definition, Wi(h) H B C Wi(h). 
Part (b) follows from proposition 12.21 (b) and (c). 

For h € H and i = 1, fc, let Wi(h) = vol(Wi(h)) > 0. For each h £ H,by proposition |2.4| applied to 
a large compact domain X whose interior contains uf =1 Wj(/i), we see that u>i(h) is a differentiable function 
so that = for all i, j = 1, k. Thus the differential 1-form rj = Ya=i Wi{h)dhi is a closed 1-form 
on the open convex set H. Since H is simply connected, there exists a C 1 -smooth function E(h) : H — > R 
so that §§- = Wi(h). 

Lemma 4.2 The Hessian matrix Hess(E) of E is positive definite for each h £ H. In particular, E is 
strictly convex and SjE : H —¥ R k is a smooth embedding. 

Proof By the same proof as in §2.3, we have for i ^ j, dwi{h)/dhj = —^pjArea(F) < if Wi{h) 
and Wj(h) share a codimension- 1 face F and it is zero otherwise. Furthermore, for each j = l,...,k, 
YLi=i dwi/dhj = d(SDi=^ w 'W) > o if Wj(h) and one of W^ + k{h) share a codimension- 1 face. It is zero 
otherwise. This shows the Hessian matrix Hess(E) = [ay] is diagonally dominated so that ay < for 
all i ^ j and an > Ylj^i \ a ij\- Thus Hess(E) is positive semi-definite. To show that it has no kernel, 
we proceed with the same argument as in the proof of corollary 13.11 The same argument shows that if 
b = (&i, bk) is a null vector for [ay], then b\ = 62 = ••• = On the other hand, there is an index % so 
that Wi(h) and one of Wj + k{h) share a codimension- 1 face, i.e, an > Yl%=i \ a ij\- Using Ylj=i a ijb\ = 0, 
we see that b\ = 0, i.e., 6 = 0. This establishes the lemma. 

Now we prove theorem [L4l as follows. Let A = {(Ai, ...,A k )\Ai > 0} and $ = \/E : H — s- A be 
the gradient map. By lemma 1431 <3? is an injective local diffeomorphism from H to A. In particular, due 
to dim(H) = dim(A), is open in A. To finish the proof that &(H) = A, since A is connected, it 

suffices to prove that &(H) is closed in A. Take a sequence of points in H so that <J>(foW) converges to 
a point a E A. We claim that a G &(Fl). After taking a subsequence, we may assume that ftW converges 
to a point in [— oo, oo] k . We first show that {h^'} is a bounded set in R k . If otherwise, there are three 

(i) 

possibilities: (a) there is j so that h - — > — oo as i — > oo, (b) there are two indices j\ and j2 so that 
lim^oo = oo and {h^} is bounded, and (c) for all indices j, limi_>oo = oo. In the first case 
(a), due to pj 6 int{conv{v\, ...,v m )) and hj is very negative, x ■ p j + /ii- < max{x • Uj/ + gj>\j' = 
1, m} for j large for all x. This implies for i large Wj(h^) = which contradicts the assumption that 
limi$(/iW) = a E i. In the case (b) that {h^} is bounded, then the sets Wj 2 (h^) lies in a compact set B. 
For i large, x ■ pj 1 + > max{x • Vj + gj,x ■ pj 2 + h^)} for all x € B. This implies that Wj 2 (h&) = 
for large i which contradicts the assumption that limj <&(/tW) = a € A. In the last case (c), since for each 
lim^oo = oo, for any compact set B, there is an index i so that B C U* =1 W^(/t^). This implies 
that the sum of the volumes Y%=i vol(W^(h^) tends to infinity which again contradicts the assumption 
limi ^(/tW) = a € A. 

Now that /iW is convergent to a point in R k , by the continuity of the map sending h to (wi(h), ...,Wk(h)) 
on R k , we see that $(h) = a. This shows h S H and a € ®(H), i.e., is closed in A. 

Hence, given any (A\, A^) E A, there exists a unique h € H so that = (^4i, ...,^4fc). Let 

u = max{x -pi + hi\i = 1, /c + m} be the PL convex function on R n and w be its dual. By corollary 2.1, 
we conclude that the vertices of w are exactly {vi,pj\i,j} with w(vi) = gi and w(pj) = —hj so that the 
discrete Hessian of w at pi, which is Wi{h) = A{. Furthermore, by proposition 12.21 the associated convex 
subdivision of w on U has exactly the vertex set {v\, Vf,,pi, ...,pk). 
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